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under  no  load  is  reported.  The  nonlinear  effects  of  magnetic  saturation  were  considered. 

The  solution  algorithm  is  a modified  version  of  an  algorithm  previously  used  in  the  analysis  of 
nonlinear,  homopolar-inductor  alternators.  The  details  of  the  algorithm  and  the  computer  programs 
implemented  on  a minicomputer  were  described  in  earlier  MERADCOM  reports. 

This  report  contains  teletypewriter  terminal  printouts  of  actual  minicomputer  runs.  A solution 
of  the  magnetic  flux  analysis  was  reached,  but  the  number  of  iterations  required  to  reach  the  solu- 
tion indicates  that  problems  remain  in  optimizing  the  numerical  analysis  used  in  the  algorithm. 
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SUMMARY 


A two-dimensional  analysis  of  the  cross  section  of  a homopolar-inductor  alter- 
nator supplying  a rectifier  load  was  developed  in  previous  investigations.  In  those 
investigations,  the  nonlinear  effects  of  saturation  of  magnetic  materials  and  nonlinear 
load  currents  were  considered.  Some  of  the  partial-difference  equations  were  para- 
bolic. The  analysis  was  applied  for  eddy  currents,  voltage  waveshapes,  and  field  excita- 
tions. An  analysis  was  made  on  a 95-kVA  homopolar-inductor  alternator,  1 15/200-V; 
3-phase;  3400-Hz;  40,800-rpm. 

The  computer  programs  which  were  developed  were  closely  associated  with  the 
decision  to  use  the  minicomputer.  Systems  of  equations  normally  require  a large 
amount  of  computer  time  and  are  programs  which  are  expensive  to  run.  Both 
algorithm  and  computer  implementation  required  investigations  having  many  runs; 
therefore,  application  of  dedicated  minicomputer  versus  large  frame  became  an 
economical  decision.  Fortran  programs  developed  for  this  investigation  are  hardware 
dependent  particularly  for  a virtual-memory  technique  for  the  storage  of  matrix  data 
in  mass-storage  devices  instead  of  a large-core  memory. 

The  present  investigation  was  limited  to  modification  of  these  previously  de- 
veloped Fortran  programs  for  the  two-dimensional  analysis  of  the  salient-pole,  synchro- 
nous alternator.  Modification  of  the  complex  interrelated  Fortran  program  and  sub- 
routines requires  detailed  and  thorough  study  and  evaluation  for  each  modification. 

This  investigation  is  limited  to  no-load  analysis  of  a model  simulating  a design  of  a 
250-kVA;  120/240-V,  400-Hz,  8000-rpm  salient-pole  aircraft  alternator.  The  dis- 
cussion of  the  analysis  emphasizes  the  interactive  characteristics  of  the  system  of  mini- 
computer programs  applied  in  the  solution  of  this  different  machine  configuration. 
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PREFACE 


This  investigation^  authorized  under  Project  1G762708AH67,  was  performed  in 
the  Electrical  Equipment  Division,  Electrical  Power  Laboratory  (formerly  Laboratory 
3000),  US  Army  Mobility  Equipment  Research  and  Development  Command.  The 
work  was  performed  under  the  direction  of  Dr.  A.  L.  Jokl,  Chief,  Electrical  Equipment 
Division.  Dr.  L.  I.  Amstutz  assisted  in  the  understanding  of  the  computer  implementa- 
tion of  the  algorithm. 
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SYNCHRONOUS  ALTERNATOR  ANALYSIS 


I.  BACKGROUND 

This  analysis  of  the  salient-pole,  synchronous  alternator  is  related  to  a similar 
analysis  performed  on  the  homopolar-inductor  alternator. 

As  in  the  homopolar-inductor-altemator  analysis,  the  numerical  method  for 
solving  the  partial-difference  equations  is  a strongly  implicit  iterative  procedure  taking 
into  account  the  periodicity  conditions  of  the  salient-pole,  synchronous  alternator. 

II.  PURPOSE 

The  investigation  demonstrates  the  analysis  of  salient-pole,  synchronous 
machinery  applying  algorithms  previously  developed  for  the  analysis  of  the  homopolar- 
inductor  alternator.  This  investigation  required  establishing  a suitable  rectangular  grid 
system  and  selecting  appropriate  boundary  conditions  for  the  specification  of  the  new- 
machine  configuration. 

HI.  PROGRAM  PRINTOUTS 

Several  problem  configurations  for  simulation  of  the  salient-pole,  synchronous 
alternator  were  tried.  Such  configurations  were  based  on  experimentation  related  to 
opposing  characteristics  of  the  simulation  and  the  solution  algorithm.  Under  load  con- 
ditions, symmetrical  boundary  conditions  should  be  specified  for  a pole  pair.  This 
configuration  provided  a grid  space  which  was  wide  compared  to  its  height.  There  is 
evidence  that  the  solution  algorithm  converges  more  slowly  when  the  grid  space  is 
wider  than  it  is  tall. 

The  result  was  to  change  to  a single-pole,  simulation  model.  The  single-pole 
models  were  considered  primarily  to  increase  the  rate  of  convergence  of  the  solution 
and  to  reduce  the  execution  time  by  specifying  a fewer  number  of  grid  points. 

IV.  PROGRAM  MODIFICATION 

The  analysis  of  homopolar-inductor  alternators  was  developed  as  a general  solu- 
tion algorithm.1  1 It  was  anticipated  that  programs  and  subroutines  previously 


A.  L.  Jokl  and  L.  I.  Amstutz,  Report  2113,”  Nonlinear  Homopolar  Inductor  Alternatonon  Nonlinear  Loads  - 
Part  I,”  USAMERDC,  Fort  Belwir,  VA.  October  1974. 
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A.  L.  Jokl  and  L.  1.  Amatutz,  Report  2114,”  Nonlinear  Homopolar  Inductor  Alternatonon  Nonlinear  Loads  - 
Part  II  ” USAMERDC,  Fort  Balvoir.  VA,  October  1974. 
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developed  for  the  homopolar-inductor  alternator  could  be  modified  to  provide  a 
solution  algorithm  for  the  salient-pole,  machine  configuration. 

Programs  and  subroutines  implemented  for  the  solution  algorithm  combine  to 
form  a complex  programming  system.  Modification  of  such  a complex  programming 
system  is  not  a trivial  exercise.  Prior  to  modification,  a time-consuming,  detailed  study 
is  required  on  each  portion  of  the  program.  Extensive  use  of  common  variables, 
multiple  use  of  certain  array  or  matrix  variables,  and  logical  branching  are  the  pro- 
gramming factors  requiring  detailed  analysis  with  modification. 

The  primary  reason  to  modify  the  programs  and  subroutines  was  to  make  the 
change-of-boundary  conditions  resulting  from  different  machine  configurations.  The 
boundary  conditions  are  applied  to  values  of  the  magnetic  vector  potential,  A.  For  the 
left  and  right  boundaries,  magnetic  vector  potentials  are  identical  for  the  homopolar- 
inductor  alternator  but  are  opposite  sign  for  the  salient-pole  alternator.  Similarly, 
lower  and  upper  boundaries  differ  to  permit  a through  flux  for  the  homopolar  con- 
figuration and  to  permit  no  through  flux  for  lower  and  upper  boundaries.  Boundary 
changes  in  the  algorithm  were  provided  by  introducing  two  new  variables: 

CS1GN  = -1.0D00 
CSYN  = 0.0D00 

Assignment  to  alternate  values  of  these  variables  and  recompiling  should  make  solu- 
tion of  the  homopolar  configuration  possible: 

CSIGN  * 1.0D00 
CSYN  = 1.0D00. 

Modification  for  the  purpose  of  boundary-condition  changes  represents  planned 
changes.  Other  planned  changes  included  modification  of  the  subroutine  SETUP  used 
by  program  UPDATE  to  change  values  of  certain  key  common  variables.  Other 
planned  changes  in  UPDATE  reduced  certain  duplicated  teletypewriter  input  used  in 
the  description  of  the  configuration.  Additional,  planned  changes  made  to  program 
FLUX  speeded  the  graphic  output  of  the  machine  configuration  on  the  CRT  terminal. 

Certain  other  modifications  of  programs  and  subroutines  provide  an  effort  to 
shorten  programs  and  to  reduce  program  execution  time.  Particular  attention  was 
given  to  program  SOLVE  and  its  functions  and  subroutines  - Subroutines  SAEB  and 
TCTB  - which  determine  reluctivity  of  iron  from  magnetic  induction  by  a table-look- 
up procedure.  Table  variables  were  changed  from  Fortran  assignment  statements  to 
data  statements.  Saving  in  execution  time  by  this  modification  was  minimal  because 
these  subroutines  were  executed  only  when  recalculation  of  reluctivities  was  required 
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instead  of  every  iteration  of  SOLVE.  Other  attempts  were  made  to  reduce  the  number 
of  calls  of  special  functions  that  require  input  and  output  from  disc  files. 

Program  modification  and  program  execution  both  produce  a large  amount  of 
teletypewriter  or  line-printer  output.  In  order  to  provide  for  a chronological  order  to 
modifications  and  program  runs,  a system  Fortran  compatible  subroutine  called  DATE 
was  devised.  A patch  was  introduced  in  the  system  program  EDITOR  at  the  beginning 
of  a file  edit.  These  system  applications  and  minor  modifications  provide  a valuable 
documentation  aid  used  on  these  programs  and  others  in  support  of  research  and  de- 
velopment projects. 

The  programs  and  subroutines  modified  during  this  investigation  are  shown  in  the 
Appendix. 

V.  DISCUSSION 

This  analysis  essentially  represents  a numerical  analysis  for  the  solution  of  partial 
differential  equations.  The  analysis  consequently  contains  those  mysterious  tricks 
associated  with  solutions  of  partial  differential  equations  and  another  set  associated 
with  numerical  analysis.  The  analysis  of  equipment  retains  certain  experimental 
characteristics  in  spite  of  all  the  efforts  to  refine  the  procedures. 

The  numerical  solution  is  provided  by  an  iterative  procedure  designed  to  give  an 
approximation  for  the  solution  within  a specified  limit  of  error.  The  procedure  does 
provide  a solution,  but  the  experimental  investigations  attempt  to  achieve  the  solution 
within  a fewer  number  of  iterations.  A sluggish  convergence  to  the  solution  is  a signifi- 
cant deterrent  regardless  of  the  size  and  speed  of  the  computer  system  which  is  applied 
to  its  solution. 

Another  problem  in  the  procedure  requires  more  emphasis  to  be  placed  upon  the 
problem  grid  than  on  easy  and  rapid  specification  of  machine  component  sizes  from  a 
drawing.  Certainly,  more  effort  is  required  to  reduce  the  tedious  and  time-consuming 
effort  in  establishing  or  significantly  modifying  the  equiyalent  model  of  the  procedure. 
Another  aspect  of  this  problem  is  that  model  definition  in  terms  of  polar  coordinates 
would  be  preferred  to  the  present  rectangular  coordinates.  The  significance  of  modi- 
fication of  the  algorithm  coordinate  system  is  an  investigation  remaining  to  be 
performed. 

The  algorithm  for  this  numerical  analysis  of  electric  machinery  is  implemented 
with  a series  of  computer  programs  supported  by  a large  number  of  functions  and  sub- 
routines. The  organization  of  programs  is  desirable  in  that  progress  of  solution  can  be 
observed  in  intermediate  results.  The  organization  provides  considerable  amount  of 


error  messages  to  guard  against  faulty  data  in  parameters  passed  in  functions  or  sub- 
routines. Disadvantages  of  the  modular  organization  include  a difficult  procedure  for 
modifying  common  variables  in  each  program  and  subroutine,  and  certain  boundary 
conditions  are  imbedded  in  subroutine  programs  and  add  to  the  complexity  of  de- 
veloping a common  solution  algorithm  when  machine  configurations  require  modifica- 
tion of  boundary  conditions. 

VI.  CONCLUSION 

It  is  concluded  that  algorithms  implemented  in  Fortran  programs  and  subroutines 
for  the  numerical  solution  of  a homopolar-inductor  alternator  provide  a solution  of  a 
salient-pole  alternator  after  correction  for  different  boundary  conditions. 
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UPDATE  AND  ORDER 


The  solution  algorithm  is  a set  of  interactive  Fortran  programs  and  subroutines. 
In  addition  to  a numerical  analysis  of  a set  of  partial  differential  equations,  this 
algorithm  is  a data-structure  problem  with  programs  for  input,  output,  and  inspection 
or  modification  of  the  data. 

The  printouts  for  several  of  the  main  programs  of  the  analyses  provide  a record  of 
the  solution.  The  interactive  nature  of  the  program  is  emphasized,  and  the  printouts 
supply  examination  of  the  machine  configuration  of  the  data  structure  and  the  status 
of  progress  toward  the  desired  solution. 

The  program  UPDATE  is  the  main  program  to  establish  the  data  structure  of  a 
new  configuration  and  to  modify  the  structure.  The  printout  for  program  UPDATE 
shows  a modification  procedure  for  any  of  13  common  variables.  Prior  to  modifica- 
tion of  subroutine  SETUP,  assignment  of  values  to  these  common  variables  required 
loading  and  execution  of  a short  Fortran  program  called  PATCH. 

This  printout  shows  the  use  of  the  Fortran  subroutine  DATE  which  marks  the 
date  on  the  printout  for  modification  of  the  data  structure. 

The  symbol  $ is  teletypewriter  prompting  for  interactive  keyboard  input.  A zero 
or  carriage  return  parameter  causes  another  printout  of  the  consolidated  variable 
listing.  A negative  parameter  causes  an  exit  to  main  program  from  SETUP  where 
options  to  change  data  at  designated  nodes  of  the  problem  grid  can  be  exercised. 

The  program  ORDER  performs  its  task  on  the  data  structure  without  any  opera- 
tor input  except  for  system  commands  associated  with  load  and  go. 
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The  printout  of  program  LIST  consists  of  two  parts.  A modification  was  made  to 
provide  a brief  heading  to  identify  the  problem,  the  date,  and  the  grid  size.  The  first 
part  is  a grid  map  that  identifies  nodes  with  grid  space  configuration.  The  first  row  of 
data  identifies  the  width  coordinate  I (1  through  IMAX).  The  first  column  of  this  part 
of  the  printout  identifies  the  height  coordinate  V (VMAX  through  1).  This  repre- 
sents data  of  the  ID  matrix  printed  on  three  pages  by  the  subroutine  IDOUT. 


The  second  printout  of  LIST  identifies  nodes  with  nine  elements  associated  with 
every  node.  This  output  comes  from  ITYPE  matrix  and  results  from  execution  of  sub- 
routine TPOUT.  Columns  2 through  10  represent  elements  of  the  nodes. 
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MAP 


Printouts  for  program  MAP  and  several  MAP  subroutines  provide  six  separate 
maps: 

1 . Boundary  Map 

2.  Conductivity  Map 

3.  Rotor  Movement  Map 

4.  Iron  Map 

5.  Current  Map 

6.  Dimensions  Check  Map 

The  negative  problem  number  in  the  heading  of  the  map  is  an  indicator  that  a 
subsequent  program  START  needs  to  be  executed.  It  is  noted  that  all  six  maps  are 
made  even  though  no  useful  information  is  provided  in  two  maps.  These  maps  provide 
the  major  comparison  of  the  data  structure  and  the  simulation  model  specifications. 

The  boundary  map  depicts  external  boundaries  with  an  E and  internal  boundaries 
or  material  changes  with  a B.  The  model  shown,  using  air  space  above  the  stator  and 
below  the  rotor,  was  intended  to  provide  zero  flux  through  stator-air  interface  above 
and  zero  flux  through  rotor-air  interface  below.  This  was  an  effort  to  provide  the  hori- 
zontal boundaries  with  configuration  of  the  model  rather  than  changing  to  software 
boundary  conditions.  The  primary  disadvantage  is  that  grid  locations  in  these  air 
spaces  require  execution  time  and,  perhaps,  solution  stability  problems  with  larger 
numbers  of  grid  points. 

The  conductivity  map  and  the  rotor  movement  map  contain  no  useful  informa- 
tion in  this  example. 

In  the  iron  map,  S indicates  SAE  4340  iron  in  the  rotor  and  T indicates  Trancor  T 
steel  in  the  stator.  A blank  indicates  air. 

The  field  current  in  the  rotor  coil  is  indicated  with  A and  X on  the  current  map. 
The  machine  was  investigated  for  no  load  so  no  currents  are  shown  for  stator  currents. 
The  stars,  or  asterisks,  provide  information  related  to  summation  of  currents  at  nodes 
adjacent  to  the  coil. 

The  dimensions  check  map  checks  for  consistent  node  spacing  between  the  grid 
nodes.  If  all  spaces  are  consistent,  the  entire  map  is  blank  except  for  coordinate 
symbols  at  head,  left  ride,  and  bottom  of  map.  An  asterisk  is  an  indication  for  a 
dimension  error.  An  error  in  this  example  occurred  in  rows  1 and  29  because  distance 
to  node  above,  H2,  and  distance  to  node  below,  H4,  were  unequal  values. 
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START  AND  SOLVE 


The  program  START  takes  the  configuration  specification  which  has  been  input 
and  checked  by  previous  programs  and  provides  an  initial  solution  to  magnetic  vector 
potentials  in  the  A matrix.  After  loading,  an  input  on  the  keyboard  of  the  negative 
value  of  problem  number  performs  the  initialization  calculations.  The  problem 
number  is  changed  to  a positive  value  as  the  indication  that  this  program  was  executed. 

The  next  output  comes  from  execution  of  the  program  SOLVE.  This  output 
indicates  that  the  conditions  for  solution  were  reached  after  iteration  149. 


The  conditions  of  solution  were  achieved,  but  the  number  of  iterations  is  wrong. 
A previous  execution  of  SOLVE  went  through  an  estimated  100  iterations  prior  to  this 
run.  This  procedure  provides  an  analysis  much  closer  to  the  acceptable  solution  than 
that  provided  by  Program  START.  The  reloading  and  restarting  of  SOLVE  execution 
were  required  to  retain  table  of  Beta  values  which  depends  on  logic  path  where  the 
value  of  Flux  is  equal  to  the  maximum  and  minimum  value  of  the  vector  magnetic 
potentials.  The  Beta  values  of  SO  are  used  as  keys  to  subroutines  to  recalculate  the 
reluctivities. 

Printout  of  asterisk  and  period  provides  location  pointers  at  the  beginning  of  the 
program  SOLVE.  During  the  interval  between  the  asterisk  and  period,  the  algorithm 
provides  a real  variable  printout  of  Beta.  If  the  Beta  printout  is  zero,  a recalculation  of 
reluctivities  is  not  done.  If  Beta  is  SO,  then  a recalculation  of  reluctivities  is  done  and 
the  parameters  of  subroutine  RELSET  are  printed  out.  It  was  noted  later  that  this 
parameter  listing  duplicates  information  shown  in  the  S-line  printout  for  each  iteration. 


Subroutine  RELSET  parameters  are: 

PREV,  PSREV,  IIDR,  IISDR,  NNDR,  NNSDR,  MXI,  MXJ. 
The  S-line  printout  for  each  iteration  is  described  in  the  chart: 

Line  1 


1. 

ITER 

Iteration  Number 

2. 

RMAX 

Maximum  Residual 
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( |,28) 

(35,25)  (18,24) 

(35,25)  l.E-04 

I .E-03  6.E-05  l.E-05  3.E-03  2.E-04  7.E-05  24R28  1 '3  0 0.E  00  3 


0.00 

28.35  24.55  I 33  12  87  18  24 

16  l.E  02  2.E  01  7.E  00  87NI2 

(36,10)  35V  I 

( 2,25)  (18,24) 

( 2,25)  l.E-04 

3.E-04  2.E-05  6.E-06  5.E-04  9.E-05  7.E-05  24R28  0 0 50  0.E 


50.00 

28.31  24.28  1 56  13  86  18  24 


86NI3 
36  V I 
(18,24) 

( 4,25)  l.E-04  „ 

8. E-03  5.E-04  l.E-04  9. E-03  l.E-03  4.E-04  24R2S  01  0 0.E  00  6 


17  3.E  02  4.E  01  l.E  01 

(19.23) 

(19.24) 


28 


0.00 


18  I.E  02  S.E  01  I.E  01  86N13 

( 1.10)  56V  I 

(22,24)  ( 0,  0) 

(22.24)  l.E-04 

4. E-03  3.E-04  6.E-05  6.E-03  6.E-04  2.E-04  24R28  02  0 0.E  00  6 

0.00 

• 

19  9.E  01  I.E  01  6.E  00  86N13 

(35.27)  56V  I 

( 1,25)  ( 0,  0) 

(35.25)  l.E-04 

5. E-04  2.E-05  6.E-06  1 .E-03  2.E-04  2.E-04  24R28  0 3 . 0 0.E  00  6 

* 

0.00 

0 

20  2.E  02  2.E  01  7.E  00  86N13 

(36.27)  56V  1 

( ?,25)  < 0,  0) 

( 2,25)  l.E-04 

2.E-04  2.E-05  6.E-06  3.E-04  8.E-05  6.E-05  24R28 
ITERATION  LI0IT.  DUMP  K RELOAD  TO  CONTINUE 
BUFFER  CLEARED.  SETUP  RECORDED . 

DATA  WITTEN  ON  TAPE 
THIS  IP  THE  DISC  VERSION  D 18JUNE73 

FORTRAN  STOP 

S//ABSD, 15, 120,65140,1000 
S//REV,  14 
S//REH,  15 
S//ABSL, I 4 
$//GO 

THIS  IS  THE  DISC  "ERPION  D 18JUNE73 

RESTART 

ENTER  ITERATION  LIMIT 
$050 

PRINT  EVERY  STEP 
$000 


1 3 0 0.E  00  6 


ITER  RMAX 

* 


RR1S  RCCR  D1AX  D?MS  DCCC  REL  OK  BA  B UB.E.C.  PT 


0.00 

26.55 


20.(3  I 40  16  83  19  24 


21  I.E  02  2.E  01  6.E  00  83NI6 

(36,29)  40V  I 

(34,24)  (19,24) 

( 2,25)  l.E-04 

9.E-05  I.E-05  4.E-06  2.E-04  3.r-05  2.E-05  2HR26 


0 0 50  0.E  00  7 


29 


50.. 90 


26.77  20.13 

*2?  2.E  P2  2.F  01 

(19.23)  ' 

(19.24) 

6.E-03  4.E-04 


2 41  17  82  19 

9.E  00 


24 


(19,24)  9.P-05 
6 .E-05  8.E-03  S.E 


82  Ml  7 
41V  2 
(19,24) 


•04  2.E-04  2(*»26.  0 1 0 0 


0.00 

*23  I .E  02  1 .E  01 
( 1,27) 

(18,25) 

2.E-03  9. E-05 


6.E  00 


(18,25)  9.E 
1 .E-05  3 .2-03  2.E 


82  Ml  7 
41V  2 
( 0,  0) 
■05 

•04  7. E-05  20R26 


0 0, 


P . PP 
IS.26 


1 1 .3$ 


24  1 .E  02  2.E  01 
(38,27) 

(17,25) 

5.E-P4  2.E-05 


2 30  18  83  19 

S.E  00 


(17,25)  9.E- 
6.E-06  6.E-04  5.E" 


22 

83N1  6 
30V  2 
(19,22) 
05 

05  4.E-05  I IRIS 


0 0 50  0, 


50.00 

16.04 


10.91 


25  I.E  02  2.E  01 

(19.23) 

(19.24) 


2 30  17  82  19 

7.E  00 


22 


£-05 


82  Ml  7 
30V  2 
(19,22) 


(19,24)  8.E- 

3.E-03  2.E-04  3.E-05  3.E-03  3.E-04  I .E-04  1 0R 1 6 0 1 


0 0, 


0.00 

26  2.p  02  2.E  0| 
( 1,27) 

(19,22) 

l.E-03  7.E-05 


6.E  00 


(19,22)  8.F 
I.E-05  2.E-03  2.E 


82N1  7 
30V  2 
( 0,  0) 
•05 

■04  I.E-04  I0R16 


0 2 


0 0, 


0.00 

2’7  S.E  02  2.E  01 
(19,20) 

(18,25) 


7.E  00 


82  Ml  7 
30V  2. 

( 0,  0) 


(18,25)  8.E-05 

l.E-03  4.E-05  9.E-0S  l.E-03  7.E-05  3.E-05  I0R1S  13  00 


E 00  7 


E 00  7 


E 00  8 


E 00  8 


E 00  8 


E 00  8 


30 


* 


0.00 


28  I .E  02 
(36,29) 
( 2,25) 


I .E  01  5.E  0M 


82  J1  7 

3 0V  2 
( 0,  0) 

(20,18)  8.E-05 

*.E-05  5.E-06  3.E-06  3. E-04  l.r-04  1 ,E-04  10R16  0 3 


0 0.E  00  8 


0.00 

*29  2.E  02 
(36,27) 
(25,24) 

5.E-05 


2. E 01  5. E 00  82N1 7 

30V  2 
( 0,  0) 

(17,17)  g.F-05 

7.E-06  3.E-06  2.E-04  7.E-05  6.E-05  10R16  1 3 0 0.E  00  8 


0.O0 

*30  1 .E  02 
( 1,27) 
( I,  9) 

4.E-05 


2 .E  01  5.E  00  82NI 7 

30'-'  2 

( 0,  0) 

( 1,  9)  8.E-03 

6.E-06  3.E-06  I.E-04  3.E-05  2.r-05  10RI6  0 0 50  O.E  00  9 


50.00 

10.94 

*31  I . E 02 
( 1,29) 
(19,24) 


7.81  1 25  17  82  18  23 

1 . E 01  6.E  00 


82*11  7 
25V  1 
(18,23) 


(33,25)  7.E-05 

9.E-04  7.E-05  2.E-05  2.E-03  3. E-04  8.E-05  7R10  01  0 0.E  00  9 


0.00 

*32  l.E  02 
(36,26) 
(36,25) 

4. E-04 


1 . E  01  5.E  00  82 vj  7 

25V  1 

( 0,  0) 

(36,25)  7.E-05 

2. E-05  4.E-06  7.E-04  6.E-05  4.E-05  7R10  12  0 0.E  00  9 


0.00 

*33  l.E  02 
(36,27) 
( 2,25) 


E-04 


1 .E  01  5.E  00  82NI 7 

25V  I 

( 0,  0) 

( 2,25)  7.E-05 

7.E-06  3.E-06  3. E-04  4.E-05  3.E-05  7RI0  0 0 50  0.E  00  10 


31 


* 


50.00 

14.73  10.71  2 25  20  79  1 9 24 

*34  l.E  02  2.E  01  6.E  00  79N20 

( 1,27)  25V  2 

( 5,25)  (19,24) 

( 5,25)  6.E-05 

5.E-P4  6.E-05  2.E-05  3.E-03  3.E-04  l.E-04  1 0R 1 4 0 1 0 0.E  00  10 

* ‘ 

0.00 

*35  l.E  02  2.E  01  5.E  00  79N20 

( 1,27)  25V  2 

( 1,25)  ( 0.  0) 

( 1,25)  6.E-05 

4.E-04  3.E-05  7.E-06  9.E-04  l.E-04  l.E-04  I 0H  I 4 0 2 0 0.E  00  10 

4 

0.00 

*36  2.E  02  2.E  01  5.E  00  79^0 

(36,27)  25V  2 

(34.24)  ( 0,  0) 

( 2,25)  6.E-05 

l.E-04  l.E-05  5.E-06  3.E-04  3.E-05  2.E-05  10R14  13  0 0.E  00  10 

* 

0.00 

*37  l.E  02  l.E  01  4.E  00  79N20 

(36,29)  25V  2 

(36.25)  . ( 0,  0) 

(15,23)  6.E-05 

4. E-05  4.E-06  2.E-06  3.E-04  l.E-04  l.E-04  10RI4  03  0 0.E  00  10 

* 

0.00 

*38  2.E  02  2.E  01  5.E  00  79N20 

( 1,27)  25V  2 

( 1,26)  ( 0,  0) 

( 1,26)  6.E-05 

5. E-05  6.E-06  3.E-06  2.E-04  6.E-05  5.E-05  I0RI4  13  0 0.E  00  10 


0.00 

*39  l.E  02  l.E  01  4.E 
( 1,27) 

( 1,24) 


(35,26)  6.E-05 


79N2  0 
25V  2 
( 0,  0) 


4.E-05  5.E-06  2.E-06  l.E-04  2.E-05  l.E-05  10114  0 0 50  0.E  00  II 


32 


50.00 


I 


25.45  15.63 


1 22  20  79  5 25 


*0  9.E  01  l.E  01  4,r  0fl 
(36,29) 

(20,25) 

( 6,25)  3.E-05 
3.E-04  3.E-05  Q.E-06  1 .E-03  2.E-kJ4  5.E-C5  15R25 


79.-J2  0 
92  V I 
( 5,25) 


1 I 0 0.  E 00  II 


0.00 


41  9.E  01  l.E  01  4 . E 00  79^0 

(36,27)  22 V 1 

( 2,25)  ( 0 p) 

( 2,25)  5.E-05 

l.E-04  7.E-06  3.E-0S  2.E-04  4.E-05  3.E-05  15R25 


0 0 50  0.E  00  12 


50.00 

13.39 


9.15  1 13  27  72  2 24 


42  l.E  02  l.E  01  4.E  00 
( 1,27) 

(20,25) 

( 3,25)  4.E-05 

2.E-04  l.E-05  5.E-06  3.E-04  4.E-05  3.E-05 


72N2  7 
13V  1 
( 2,24) 


9RI3  1 1 0 0.E  00  12 


0.00 


43  9.E  01  l.E  01  4.E  00  72  '12  7 

( 1,27)  13V  i 

(13,18)  ( 0,  0) 

(35,25)  4.E-05 

4.E-05  4.E-06  2.E-06  l.E-04  5.E-05  4.E-05  9R13  0 0 50  0.E  00  13 

BUFFER  CLEARED.  SETUP  RECORDED. 

DATA  WRITTEN  ON  TAPE 

THIS  IS  THE  DISC  VERSION  D 18JUNE73 


FORTRAN  STOP 

S//ABSD, 15, 120,65140,1000 
PRESETTING  SYSTEM 
S//RESET 

RESETTING  SYSTEM 
S//ATT, 4 , TTY 

P//ATT, 14, LT7, RELOAD, DISC  25  FEB  76 
S//ATT, 15, LT8, RELOAD, DISC 
S//ABSL, 14 
S//GO 

THIS  IS  THE  DISC  -’ERSION  D I8JUNE73 
RESTART 
$000 

ENTER  ITERATION  LIMIT 
$100 

PRINT  EVERY  STEP 


33 





i te^  r^m  Rr>r  'ney  occc  «» f;l  ok  ? • ? up.e.c.  pt 

* 


5 P.PP 
P.26 


.67  | ||  ??.  71  14 


4-*  1 .E  O'?  | .2  pi l 4.E  PP  7IN23 

(36.27)  ||W  1 

( 2 ,25)  ( I a ,25) 

( 2,2 5)  4.S--05 

9.E-P5  9.E-06  4.E-06  3.E-IJ4  8.E-P5  7.E-P5  PR  3 o | 0 P.E  Bf!  13 

k 

PI. PIP 

45  I .E  0?  1 .E  PI  I 4.E  PIP  7TI28 

(36.28)  I IV  I 

(14,17)  ( P,  P) 

(36,23)  4 ,E**>5 

6. E-P5  6.E-P6  3.E-06  I .E-04  2.E-05  I .E-P5  PR  8 1 2 0 P.E  PP  13 

i 

P.PP 

*46  7.E  PI  8.E  PP  3.E  PP  71N28 

( 1,29)  1 I V 1 

( 1,25)  ( 0,  P) 

(27,  9)  4.E-05 

2. E-P5  3.E-P6  1.E-P6  2.E-04  f.E-04  9.E-05  PW  8 0 3 P P.E  00  13 

P.PP 

*47  I .E  P2  1 ,E  01  4.E  00  7PJ28 

( 1,27)  IIV  I 

( 1,24)  ( 0,  0) 

( 1,24)  4.E-05 

7. E-05  5.E-06  2.E-06  2.E-P4  4.E-05  3.E-05  PH  8 I 3 P P.E  00  13 

k 

0.PP 

*48  9.E  01  I .E  01  3.E  00  71  IB  8 

(36,27)  -IIV  I 

(36,  9)  ( 0,  0) 

(36,  9)  4.E-05 

3. E-P5  4.E-06  2.S-06  7.E-05  2.E-05  l.E-05  PR  8 0 P 50  0.E  00  14 

k 

5P.  00 

10.49  3.57  I 11  23  76  19  18 

*49  7.E  01  7.E  00  3.E  00  76N23 

(36.29)  IIV  I 

(36,25)  (19,18) 

(33,25)  3.E-05 

5.E-05  7.E-06  3.E-06  3.E-04  3.E-05  2.E-05  3R 1 0 II  0 0.E  00  14 


34 


* 


<*.00 


5fl  6.E  O! 

7.E  00 

3."  00 

76 1C  3 

( 1,2?) 

1 IV  1 

( 1,05) 

(35,  9) 

3 .c-05 

( 0,  0) 

4.E-05 

3. "-06 

2.E-06  6.1-M5 

3.1-05 

2.E-05  3R10 

0 

0 

50 

0.E 

0P 

15 

50.00 

8.48 

3.7? 

1 10  27  72 

2 24 

51  7.E  01 

8.E  00 

3.E  00 

72 1C  7 

( 1,27) 

10V  I 

( 3 ,05) 

( 3,25) 

3.E-05 

( 2,24) 

7. FI-05 

7.E-06 

3.E-06  3. '•-04 

4.E-05 

2.E-05  3R  8 

1 

1 

0 

0.E 

00 

15 

0.00 

50  6.E  01 

7.E  00 

3.E  00 

72N2  7 

(36,09) 

I0v  1 

(36,25) 

(36,25) 

3.E-05 

( 0,  0) 

4.E-05 

3.E-06 

2.E-06  1.E-W4 

4.E-05 

3.E-05  3R  8 

P 

2 

0 

0.E 

00 

15 

0.00 

53  7.E.01 

8.E  00 

3.E  00 

72  “C  7 

(36,27) 

10V  1 

(28,24) 

( 7,24) 

3.E-B5 

( 0,  0) 

3.E-05 

l 

3.E-06 

l.E-06  l.E-04 

7.E-05 

6.E-05  3R  8 

0 

3 

0 

0.E 

00 

15 

0.00 

54  8.E  01 

9.E  00 

3.E  00 

72NB7 

( 1,27) 

10V  1 

( 1,23) 

( 1,  9) 

3.E-05 

( 0,  0) 

3.E-05 

r 

3.E-06 

l.E-06  7.E-05 

I.E-05 

8.E-06  3R  8 

1 

3 

0 

0.E 

00 

15 

0.00 

55  6.E  01 

6.E  00 

2.E  00 

72N2  7 

( 1,29) 

10V  1 

( 1,29) 

(33,22) 

3.E-05 

( 0,  0) 

2.E-05 

2.E-06 

l.E-06  2.E-04 

9.E-05 

7.E-05  3R  8 

0 

3 

0 

0.E 

00 

15 

0.00 

56  9.E  01 

1 . E 01 

3.E  00 

72N2  7 

(36,27) 

10V  I 

(36,  9) 

(36,11) 

3.E-05 

( 0,  0) 

3.E-05 

4.E-06 

l.E-06  8.E-05 

3.E-05 

3.E-05  3R  8 

1 

3 

0 

0.E 

00 

15 

: ' 


35 


0.00 


57  6.E  01  9.E  00  3.E  00 
(36, 1 0) 

(36,  9) 

( 1,20)  3. E-05 

2. E-05  3.E-06  l.E-06  R.E-05  I .E-05  9.E- 


72  N2  7 
10V'  I 
( 0,  0) 

06  3R  8 T 0 50  0.E  00  16 


50.00 

11  .38  6.47  i 11  20  79  29  25 

*58  5.E  01  6.E  00  3.F.  00  79, C0 

( 1,29)  II V 1 

(29,25)  (29,25) 

(29,25)  2. E-05 

4. E-05  5.E-06  3.E-06  2.E-04  3. E-05  2. E-05  6” I I 1 1 


0 0.E  00  16 


0.00 

*59  4.E  01  6.E  00  2.E  00  79N20 

(35,29)  IIV  1 

( I,  9)  ( 0,  0) 

(35,  9)  2. E-05 

1 .E-05  2.E-06  I.E-06  5. E-05  3. E-05  2, E-05  6RI 1 0 0 50  0.E  00  17 


50.00 

6.92 


1.79  1 9 26  73  22  18 


60  5.E  01  6.E  00  3.E  00  73926 

(36,27)  9V  I 

(25,18)  (22,18) 

( 8,25)  2. E-05 

3. E-05  4.E-06  2.E-06  l.E-04  2. E-05  2. E-05  1R  6 I I 0 0.E  00  17 

» 

0.00 

*61  4.E  01  6.E  00  2.E  00  73  92  6 

(36,29)  9V  I 

(10,24)  ( 0,  0) 

( 4,24)  2. E-05 

2. E-05  2.E-06  I.E-06  6. E-05  3.F-05  3.E-05  IR  6 0 2 0 0.E  00  17 


0.00 

*62  9.E  01  7.E  00  2.E  00  75926 

( 1,27)  9 V l 

( I,  9)  ( 0,  0) 

(26,  9)  2. E-05 

2. E-05  2.E-06  I.E-06  l.E-04  7. E-05  6. E-05  IR  6 0 3 0 0.E  00  17 

fc 

0.00 

*65  6.E  01  8.E  00  2.E  00  73N26 

( 1,10)  9 V I 

( 1,24)  ( 0,  0) 

( 1,24)  2.E-05 

6, E-05  4.E-06  I.E-06  8. E-05  I.E-03  7.E-06  IR  6 I 3 0 0.E  00  17 


36 


1 


0.00 


64  3.E  01  5.E  00  2.E  00  73*6 

(36,29)  9V  I 

(36,  9)  ( Pf  o) 

( 8,21)  2.E-05 

I.E-05  2.E-06  9.E-07  l.E-04  7.E-05  6.E-05  IR  6 0 3 0 0.E  00  17 


0.00 


65  7.E  01 
(36,10) 
(36,  9) 


9.E  00  2.E  00 


73*6 
9V  1 
( 0,  0) 

(36, 10)  2.E-05 
3.E-05  3.E-06  l.E-06  9.E-05  3.E-05  3.E-05  IR  6 


3 0 0.E  00  17 


0.00 


66  5.E  01 
( 1,10) 
( 1,24) 


7.E  00  2.E  00 


73*6 
9V  1 
( 0,  0) 

(36,  9)  2.E-05 

4.E-05  3.E-06  l.E-06  6.E-05  l.E-05  7.E-06  IR  6 I 3 0 0.E  00  17 


0.00 


67  3.E  01 
(35,29) 
( 1,  9) 


5.E  00  2 ,E  00 


73*6 
9 V 1 
( 0,  0) 

(35,  9)  2.E-05 

l.E-05  2.E-06  9.E-07  4.E-05  l.E-05  l.E-05  IR  6 0 0 50  0.E  00  18 


50.00 

5.10 

68  3.E  01 
(36,29) 
(13,18) 


2.88  I 1)  19  80  25  17 

5.E  00  2.E  00 


80N19 
11  V 1 
(25,17) 


(25,18)  2.E-05 
8.E-05  7.E-06  3.E-06  i.E-04  3.E-05  2.E-05  2R  5 0 I 


0 0.E  00  18 


0.00 


69  4.E  01 
(36,27) 
(14,17) 


3.E  00  2.E  00 


80  41  9 
1 1 V 1 
( 0,  0) 

(36,13)  2.E-05 

2.E-05  3.E-06  l.E-06  4.E-05  2.E-05  l.E-05  2R  5 1 2 0 0.E  00  18 


0.00 


70  3.E  01 
< 1,27) 
( I,  9) 
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ITERATION  LniT.  DUMP  S RELO'D  TC  CONTINUE 
BUFFER  CLEARED.  SETUP  RECORDED. 

DATA  WRITTEN  ON  TAPE 

THIS  IS  THE  DISC  VERSION  D I8JUNE73 

FORTRAN  STOP 

S//ABSD, 13, 120,65140, 1000 
5//REW, 14  25  FEB  76 
$//REV, 1 5 
S//ABSL, I 4 
S//GO 

THIS  IS  THE  DISC  VERSION  D I8JUNE73 
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OTFN 


The  following  programs  are  run  at  time  during  the  execution  of  SOLVE  program. 
The  output  printout  supplies  additional  field  condition  indicators  regarding  the 
suitability  of  the  solution  algorithm. 

Program  OTFN  maps  field  conditions  by  printout  on  line  printer  or  teletypewriter 
printer.  This  analysis  was  performed  only  at  location  1 , and  problems  associated  with 
rotor  motion  and  location  2 were  not  considered.  In  addition,  eddy  currents  and  eddy- 
current  losses  were  not  considered  so  maps  related  to  these  conditions  are  not  shown. 
The  maps  are  shown  in  order  as  requested  from  the  main  program. 

The  main  program  OTFN  recalculates  residual  by  call  of  Subroutine  RES.  The 
R matrix  is  used  as  scratch  matrix  as  residual  values  of  each  node  must  be  recalculated 
before  printout.  The  value  of  R matrix  is  multiplied  by  Scale  and  changed  to  integer 
values  for  printing.  The  scale  of  1000  shown  in  heading  is  multiplier  so  values  in 
residual  map  must  be  divided  by  1000  for  true  values. 

The  clustering  groups  of  negative  residuals  in  the  left-side  region  of  air  and  posi- 
tive residuals  in  the  right-side  region  of  air  instead  of  intermixed  positive  and  negative 
values  are  phenomena  which  remain  unexplained.  This  clustering  of  negative  and  posi- 
tive residuals  in  which  maximum  residual  values  occur  either  in  the  first  or  last  columns 
is  frequently  observed  in  the  iteration  printout  of  the  program  SOLVE. 

Print  B provides  an  option  to  print  a magnetic  induction,  or  B,  map.  With  this 
printout  is  a double  conversion  from  metric  to  English  units  and  from  real  to  integer 
values.  During  this  later  conversion,  a system  error  message  is  encountered  when  the 
value  of  B exceeds  32,767,  the  maximum  integer  value  for  single-precision,  16-bit 
integer  variable.  Although  the  error  is  not  fatal,  it  does,  in  addition,  provide  a highly 
saturated  flux  condition  especially  for  no-load  conditions  simulated  on  the  alternator. 

The  print  field  angle  prints  a grid  map  of  field  angles  in  degrees.  The  arc  tangent 
function  of  the  X-component  of  the  magnetic  induction  provides  the  angle  subtended 
from  the  vertical.  A negative  angle  means  the  direction  of  the  X-component  is  to  the 
left  of  vertical,  and  a positive  angle  means  the  X-component  of  B lies  to  the  right.  The 
angle  along  iron-air  interfaces  indicates  flux  at  these  interfaces  is  at  right  angles  which 
is  the  desired  condition.  The  field  angle  at  vertical  boundaries  from  iron  is  nearly  90 
degrees  which  is  again  the  proper  field  condition  for  solution  of  the  present  alternator 
configuration. 

The  relative  permeability  data  is  stored  in  DMAT  matrix  as  the  ratio  of  reluctivity 


SI 


* 

of  space,  0.7957747E+06,  to  relative  reluctivity  of  the  material.  The  relative  per- 
meability is  obtained  by  taking  the  reciprocal  of  the  DMAT  variable,  multiplying  by 
the  reluctivity  of  air,  and  converting  from  a real  to  an  integer  variable.  The  relative 
permeability  of  the  stator  steel  is  in  the  5000  range  for  unsaturated  portions  of  the 
path.  The  value  of  permeability  reduces  for  saturated  grid  spaces.  Similarly,  the  rela- 
tive permeability  of  the  rotor  iron  is  in  the  300  range.  The  relative  permeability  of  air 
grid  spaces  has  value  of  unity. 

The  vector  potential  map  in  units  of  webers/meter  consists  of  a vector  in  Z-plane. 
The  positive  value  of  vector  potential  can  be  considered  a vector  out  of  the  page,  or 
plane,  of  alternator  configuration.  On  the  opposite  boundary,  the  negative  value  can 
be  considered  a vector  into  the  page. 

A similar  group  of  printouts  may  be  selected  at  rotor  in  location  2.  However, 
this  option  was  not  exercised  for  this  example,  and  certain  variables  requiring  setup  for 
material  motion  were  never  set.  These  printouts  serve  mainly  to  indicate  values  which 
can  be  associated  with  graphic  representations  shown  in  the  next  program. 
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614 

768 

93  8 

1 123 

1016 

645 

407 

10 

176 

275 

408 

557 

719 

892 

1075 

947 

514 

146 

9 

81 

22  7 

377 

534 

700 

873 

1053 

91  1 

440 

1 70 

8 

16292 

16311 

16349 

16408 

16492 

16608 

16792 

17014 

17200 

17764 

7 

16310 

16330 

16369 

16430 

16515 

16633 

16800 

16979 

17159 

17492 

6 

16326 

1634  7 

16387 

16448 

16534 

16649 

16795 

15948 

17100 

17241 

5 

16349 

16370 

16411 

164  74 

16559 

16664 

16781 

16900 

16985 

16911 

16364 

16385 

1642  7 

16489 

16570 

16664 

16758 

16826 

16819 

16652 

3 

50 

51 

54 

58 

64 

71 

79 

88 

97 

106 

2 

28 

30 

34 

40 

47 

55 

64 

72 

81 

89 

f 

6 

12 

21 

29 

38 

47 

. 57 

65 

74 

82 

INDUCTION  IN  GAUSS  LOCATION  = I PROBLEM  L J5  1 


J 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

29 

900 

93  7 

965 

985 

998 

1007 

1013 

1015 

1015 

1013 

28 

899 

936 

964 

985 

998 

1007 

1013 

1015 

1015 

1013 

27 

956 

984 

1002 

1013 

1020 

1024 

•026 

102  6 

1026 

1026 

26 

1 104 

1096 

1089 

1076 

1061 

1051 

104  6 

1043 

1043 

1046 

25 

17831 

17091 

16091 

14351 

12052 

10285 

9M6 

4506 

4504 

91  A3 

24 

17085 

17739 

15312 

8774 

14075 

13742 

7443 

13677 

13676 

7443 

23 

504 

17352 

18415 

317 

19201 

19311 

379 

193  79 

193  76 

3 78 

22 

655 

1822  6 

18012 

368 

19238 

19323 

408 

19394 

19391 

408 

21 

914 

18221 

18290 

644 

19226 

192  84 

63  7 

19333 

19331 

63  7 

20 

1927 

17871 

17893 

2070 

13663 

18702 

2072 

1869? 

18690 

2072 

19 

2949 

17541 

17389 

3326 

18142 

18135 

3326 

18096 

18096 

3326 

18 

4161 

17852 

19144 

2630 

18709 

18775 

2358 

1862  0 

18620 

2358 

17 

2 0933 

202  3 6 

20280 

11049 

173  78 

17737 

8318 

1 7007 

17007 

8319 

16 

19337 

19095 

17065 

16916 

16348 

15351 

15371 

15209 

15210 

15371 

15 

17604 

16921 

17094 

16676 

16364 

162  93 

I6I50 

16114 

16115 

16150 

I A 

16935 

17100 

1683  7 

16780 

16701 

16610 

16578 

16553 

16553 

16577 

13 

17041 

16963 

16954 

16896 

1684  7 

16818 

16794 

16782 

1678? 

16793 

12 

1 7054 

17054 

17017 

16987 

16958 

16935 

16922 

1691? 

1691? 

16921 

11 

17134 

17091 

17082 

17048 

17022 

17007 

16983 

16967 

16966 

16982 

10 

1 722  5 

17190 

17159 

1 7099 

17045 

16997 

1694  6 

16918 

16918 

16945 

9 

17752 

17481 

17253 

17073 

16929 

16814 

16739 

16703 

1670? 

16738 

8 

18946 

17549 

17075 

16746 

16506 

16348 

16260 

1622  1 

16220 

16258 

7 

17554 

17284 

16588 

16119 

15820 

15559 

15383 

152  84 

15284 

15381 

6 

17080 

16567 

16051 

154  73 

15035 

144  83 

14015 

13706 

13  705 

140M 

5 

16533 

15904 

15154 

14150 

12913 

11546 

10356 

9684 

9684 

10354 

4 

16206 

15325 

13890 

11938 

9665 

7306 

5106 

3565 

3564 

5103 

3 

114 

119 

121 

123 

124 

125 

126 

127 

127 

126 

2 

96 

102 

107 

1 1 1 

114 

116 

117 

118 

113 

117 

1 

89 

96 

101 

106 

109 

112 

1 13 

114 

114 

1 13 

! 


L 


1 


INDUCTION  IN  GAUSS  LOCATION  = I PPOPLE!»l  LJ51 


J 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

29 

1008 

999 

985 

966 

938 

901 

854 

795 

725 

647 

28 

1008 

999 

985 

965 

93  7 

900 

851 

791 

719 

639 

27 

1024 

1020 

1014 

1003 

984 

957 

918 

863 

795 

719 

26 

1051 

1061 

1076 

1089 

1097 

1104 

1 107 

1 06S 

975 

903 

25 

102  80 

12050 

14352. 

16091 

17092 

17833 

19073 

2003  8 

20234 

90045 

24 

13  721 

14099 

8775 

15316 

17744 

17088 

17860 

20177 

20541 

19505 

23 

192  73 

19251 

317 

18423 

17361 

504 

84  78 

109  72 

2119 

5482 

22 

19289 

19282 

368 

18018 

1823  6 

654 

10254 

10552 

2115 

2910 

21 

192  5 0 

19269 

644 

18297 

18230 

914 

11011 

10982 

21^5 

3166 

20 

18685 

18687 

2072 

17900 

17878 

1928 

1 1229 

10948 

2383 

4013 

19 

18134 

18149 

3327 

17397 

17549 

2950 

10773 

10725 

2604 

2786 

18 

18775 

18715 

2631 

19155 

17858 

4 164 

13316 

9716 

3 170 

2508 

17 

1 773  7 

173  82 

11056 

20292 

20233 

20924 

21817 

11800 

6264 

3191 

16 

15351 

16353 

16924 

17062 

19093 

19334 

2 733 

1188 

1959 

1883 

15 

16293 

16366 

16677 

17090 

16928 

17604 

1699 

1317 

1559 

1581 

14 

16609 

16701 

16779 

16837 

17105 

1693  7 

1245 

1187 

1375 

1408 

13 

16817 

16845 

16894 

16959 

16964 

17042 

930 

998 

1231 

1285 

12 

16933 

16954 

16985 

17024 

17054 

17052 

661 

81  1 

1110 

1 192 

1 1 

1 7004 

17018 

17045 

17091 

17092 

17130 

406 

645 

1015 

1 122 

10 

16994 

1 7040 

17097 

17167 

17191 

1 7220 

146 

514 

947 

1074 

9 

16811 

16924 

1 7069 

17260 

17483 

17749 

171 

440 

91  1 

1052 

8 

16345 

16502 

16742 

1 7080 

1 7552 

18943 

17758 

17194 

1 7008 

16787 

7 

15556 

15816 

16116 

1S590 

172  88 

17556 

17493 

17159 

16980 

16800 

6 

14481 

15032 

154  70 

1 5051 

16569  • 

1 7084 

17246 

17104 

16950 

16796 

5 

1 1544 

1291 1 

14148 

15152 

15901 

16530 

16908 

16982 

16898 

16780 

4 

73  03 

9661 

11934 

13886 

15322 

16204 

16649 

16816 

16824 

16756 

3 

125 

124 

123 

121 

119 

114 

106 

97 

88 

79 

2 

116 

114 

1 1 1 

107 

102 

96 

8? 

«1 

72 

64 

1 

1 12 

109 

106 

101 

9G 

89 

82 

74 

65 

57 

i 


{ 


INDUCTION  IN  GAUSS 


LOCATION  : 1 


PROBLEM  LJ5 I 


J 

31 

32 

33 

34 

35 

36 

29 

561 

471 

378 

286 

205 

145 

28 

550 

456 

359 

259 

163 

75 

27 

639 

555 

475 

405 

351 

318 

26 

845 

757 

685 

646 

608 

583 

25 

20233 

20401 

20135 

201  79 

203  82 

20140 

24 

20187 

20791 

19878 

20067 

20786 

19977 

23 

7726 

2351 

6551 

7086 

2252 

6796 

22 

3957 

2351 

2925 

3094 

2256 

292  7 

21 

3225 

2353 

2510 

24  72 

2259 

2403 

20 

3029 

2358 

2775 

2596 

22  61 

2555 

19 

852 

1180 

1038 

703 

1060 

778 

18 

1390 

1187 

983 

808 

887 

772 

17 

1813 

1301 

1028 

890 

846 

790 

16 

1525 

1228 

1018 

891 

820 

779 

15 

1340 

1 125 

958 

842 

769 

732 

14 

1204 

1022 

875 

765 

691 

654 

13 

1095 

925 

783 

671 

593 

553 

12 

1006 

838 

692 

572 

483 

434 

11 

93  7 

767 

613 

479 

371 

3 05 

10 

891 

718 

556 

407 

2 73 

174 

9 

872 

S98 

533 

3 76 

223 

76 

8 

16605 

16490 

16406 

16348 

16310 

162  92 

7 

16633 

16514 

16429 

163  69 

16330 

16310 

6 

16649 

16534 

16448 

16386 

16346 

16326 

5 

16662 

16557 

164  73 

16411 

16369 

16349 

4 

16662 

16563 

16488 

16426 

16385 

16364 

3 

71 

64 

58 

54 

51 

50 

2 

55 

47 

40 

34 

30 

28 

1 

48 

39 

30 

21 

13 

6 

60 


PRINT  FIELD  ANGLE 
$ 


FIELD 

ANGLE 

IN  DEGREES 

LOC  i 

TION  = I 

1 

PROBLEM 

LJ51 

J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

29 

71 

42 

28 

20 

15 

12 

9 

7 

6 

28 

-49 

-20 

-12 

-8 

-6 

-5 

-4 

-3 

-2 

27 

-81 

-63 

-49 

-39 

-32 

-25 

-21 

-17 

-13 

26 

-86 

-75 

-64 

-57 

-48 

-40 

-34 

-28 

-21 

25 

-88 

-89 

-87 

-89 

-88 

-86 

-88 

-87 

-84 

24 

-77 

-89 

-75 

-79 

-89 

-75 

-82 

-88 

-69 

23 

-50 

-89 

-47 

-53 

-87 

-40 

-61 

-85 

-14 

22 

-74 

-89 

-67 

-82 

-89 

-46 

-72 

-88 

-10 

21 

-88 

-89 

-80 

-83 

-89 

-53 

-53 

-84 

-9 

20 

-82 

-89 

-87 

-76 

-86 

-70 

-52 

-62 

-8 

19 

-53 

-88 

-61 

-44 

-75 

-52 

-27 

-38 

-2 

18 

-75 

-88 

-81 

-66 

-76 

-66 

-46 

-45 

-2  7 

17 

-87 

-39 

-87 

-86 

-87 

-86 

-85 

-88 

-72 

16 

-88 

-85 

-83 

-82 

-30 

-75 

-64 

-42 

-79 

15 

-86 

-82 

-77 

-73 

-69 

-63 

-54 

-52 

-81 

14 

-85 

-77 

-70 

-64 

-59 

-53 

-4  7 

-50 

-71 

13 

-83 

-72 

-63 

-55 

-49 

-44 

-39 

-43 

-62 

12 

-80 

-65 

-54 

-45 

-39 

-3^ 

-30 

-34 

-52 

1 1 

-75 

-54 

-41 

-33 

-28 

-24 

-20 

-23 

-38 

10 

-63 

-3  5 

-24 

-18 

-15 

-12 

-10 

-11 

-19 

9 

-9 

-3 

-2 

-2 

-1 

-1 

0 

0 

4 

8 

89 

89 

88 

88 

87 

87 

86 

87 

88 

7 

89 

89 

88 

88 

88 

f 7 

87 

88 

89 

6 

8Q 

39 

89 

88 

88 

88 

88 

88 

-89 

5 

89 

89 

89 

89 

88 

88 

89 

89 

89 

4 

89 

89 

89 

89 

89 

89 

89 

89 

89 

3 

87 

78 

69 

61 

54 

48 

43 

37 

32 

2 

85 

68 

54 

44 

36 

30 

26 

22 

18 

1 

62 

25 

15 

10 

8 

6 

5 

4 

3 

61 


FIELD  ANGLE  IN  DEGREES 


LOCATION  = 1 


PROBLEM  LJ5 1 


• J 

1 1 

12 

13 

14 

15 

16 

17 

18 

1R 

20 

29 

3 

2 

2 

1 

28 

-1 

-1 

27 

-7 

-5 

-4 

-2 

-1 

-1 

26 

-1 1 

-7 

-5 

-3 

-1 

-1 

25 

-83 

-83 

-78 

-73 

-73 

-63 

-34 

-17 

16 

34 

24 

-84 

-56 

-3  7 

-64 

-2  7 

-6 

-15 

-10 

10 

15 

23 

-66 

1 

2 

-4 

0 

0 

-2 

0 

0 

2 

22 

-75 

-2 

-1 

-26 

0 

0 

-5 

0 

0 

6 

2 1 

-66 

-1 

-1 

-28 

0 

0 

-5 

0 

0 

5 

20 

-32 

-1 

-4 

-11 

0 

-2 

-2 

1 

-1 

2 

1 9 

-20 

-1 

0 

-4 

18 

-47 

-10 

-4 

-13 

-2 

0 

-1 

-1 

1 

1 

1 7 

-57 

-23 

-30 

-3  6 

-4 

-19 

-14 

7 

-7 

4 

16 

-12 

-2  0 

-23 

-17 

-16 

-12 

-6 

-2 

2 

6 

15 

-8 

-9 

-8 

-8 

-6 

-4 

-3 

-1 

1 

3 

1 4 

-3 

-4 

-4 

-3 

-3 

-2 

-1 

0 

0 

1 

13 

-3 

-2 

-2 

-2 

-1 

-1 

12 

-2 

-1 

-1 

-1 

-1 

1 1 

-1 

10 

0 

9 

2 

4 

4 

4 

3 

? 

1 

0 

0 

-1 

8 

45 

23 

16 

1 1 

8 

5 

3 

1 

-1 

-3 

7 

66 

45 

30 

22 

16 

1 1 

6 

2 

-2 

-6 

6 

73 

59 

44 

34 

26 

19 

I 1 

4 

-4 

-1  1 

5 

80 

71 

60 

51 

42 

32 

21 

7 

-7 

-21 

4 

86 

83 

79 

75 

69  . 

62 

49 

21 

-21 

-49 

3 

21 

16 

13 

9 

7 

5 

2 

0 

0 

-2 

2 

12 

10 

8 

6 

4 

3 

1 

0 

0 

-1 

12  2 11 


FIELD 

ANGLE 

IN  DEGREES 

LOCATION  = 1 

PROBLEM 

LJ5 1 

J 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

29 

0 

-1 

-1 

-2 

-2 

-3 

-4 

-6 

-7 

-9 

2 8 

0 

0 

0 

0 

1 

1 

1 

2 

3 

4 

27 

1 

1 

2 

4 

5 

7 

10 

13 

17 

21 

26 

1 

1 

3 

5 

7 

1 1 

15 

21 

28 

» 

25 

63 

73 

73 

78 

83 

83 

80 

84 

87 

88 

24 

6 

27 

64 

37 

56 

8^ 

70 

69 

88 

82 

23 

0 

0 

4 

-2 

-1 

65 

12 

14 

85 

61 

22 

0 

0 

26 

1 

2 

75 

6 

10 

88 

2 

21 

0 

0 

28 

1 

1 

66 

6 

9 

84 

54 

20 

2 

0 

11 

4 

1 

32 

10 

8 

62 

S 

1? 

0 

0 

4 

0 

1 

20 

3 

2 

38 

27 

18 

0 

2 

13 

4 

10 

48 

21 

27 

45 

46 

17 

19 

4 

36 

30 

23 

57 

67 

72 

88 

85 

16 

12 

16 

17 

23 

20 

12 

58 

79 

42 

64 

15 

5 

6 

8 

8 

9 

8 

79 

81 

52 

54 

14 

2 

3 

3 

4 

4 

3 

89 

71 

50 

* 

13 

1 

1 

2 

2 

2 

3 

84 

62 

43 

39 

12 

0 

1 

1 

1 

1 

2 

80 

52 

34 

0 

1 1 

0 

0 

0 

0 

0 

1 

75 

38 

23 

20 

10 

0 

0 

0 

0 

0 

0 

62 

19 

1 1 

0 

9 

-2 

-3 

-4 

-4 

-4 

-2 

-70 

-4 

8 

-5 

-8 

-11 

-16 

-23 

-45 

-87 

- 88 

-87 

-86 

7 

-1  1 

-16 

-22 

-30 

-45 

-66 

-85 

89 

-88 

-87 

6 

-19 

-26 

-34 

-44 

-59 

-73 

-85 

-89 

-88 

-88 

5 

-32 

-42 

-51 

-60 

-71 

-80 

-86 

-89 

-39 

-88 

4 

-62. 

-69 

-75 

-79 

-83 

-86 

-88 

-89 

-89 

-89 

3 

-5 

-7 

-9 

-13 

-16 

-21 

-26 

-32 

-37 

-43 

2 

-3 

-4 

-6 

-8 

-10 

-12 

-15 

-18 

-22 

-K 

1 

0 

0 

-1 

-1 

-2 

-2 

-3 

-3 

-4 

-5 

63 


-V 


FIELD 

ANGLE 

IN  DEGREE? 

LOCATION  = 

1 PROBLEM  LJ5 1 

J 

31 

32 

33 

34 

35 

36 

29 

-12 

-15 

-20 

-27 

-41 

-68 

28 

5 

6 

8 

11 

19 

46 

27 

25 

31 

39 

49 

63 

80 

26 

40 

48 

57 

63 

74 

86 

25 

86 

88 

89 

87 

89 

88 

24 

75 

89 

79 

76 

89 

77 

23 

40 

87 

53 

47 

89 

50 

22 

46 

89 

82 

67 

89 

74 

21 

53 

89 

83 

80 

89 

88 

20 

70 

86 

76 

87 

89 

82 

19 

52 

75 

44 

61 

88 

52 

18 

66 

76 

66 

8! 

88 

74 

17 

86 

87 

86 

87 

89 

87 

16 

75 

80 

82 

83 

86 

89 

15 

63 

69 

73 

77 

82 

87 

14 

53 

59 

65 

70 

78 

85 

13 

44 

49 

56 

63 

72 

84 

12 

34 

39 

45 

54 

66 

81 

1 1 

28 

33 

41 

55 

76 

10 

12 

15 

18 

24 

36 

65 

9 

• 1 

1 

2 

2 

4 

10 

8 

-87 

-87 

-88 

-88 

-89 

-89 

7 

-87 

-88 

-88 

-88 

-89 

-89 

6 

-88 

-88 

-88 

-89 

-89 

-89 

5 

-88 

-88 

-89 

-89 

-89 

-89 

4 

-89 

89 

-89 

-89 

-89 

-89 

3 

-48 

-54 

-61 

-69 

-76 

-85 

2 

-30 

-36 

-44 

-54 

-65 

-82 

1 

-6 

-8 

-10 

-15 

-23 

-52 

64 


PRINT  PORTABILITY 
$ 

RELATIVE  PERMEABILITY  LOCATION  = 1 PROBLEM  LJ51 


J 

1 

2 

3 

4 

5 

29 

1 

1 

1 

1 

1 

28 

1 

1 

1 

1 

1 

27 

1 

1 

1 

l 

1 

26 

1 

1 

1 

1 

1 

25 

30 

24 

29 

30 

23 

24 

-3  7 

-18 

-34 

-38 

17 

23 

5687 

1 

5745 

5620 

1 

22 

3920 

1 

4056 

3912 

1 

21 

3392 

1 

3451 

3485 

1 

20 

3525 

1 

3566 

3745 

1 

1 9 

l 

1 

1 

1 

18 

1 

1 

1 

1 

1 

17 

1 

1 

1 

1 

1 

16 

l 

1 

l 

1 

15 

1 

1 

1 

1 

1 

14 

1 

1 

l 

1 

1 

13 

1 

1 

1 

1 

1 

12 

1 

1 

1 

1 

l 

1 1 

1 

1 

1 

1 

1 

10 

1 

1 

l 

l 

1 

9 

1 

1 

1 

1 

1 

8 

280 

278 

274 

269 

261 

7 

279 

276 

2 72 

266 

259 

6 

277 

275 

271 

265 

257 

5 

2 75 

272 

268 

262 

255 

4 

273 

271 

267 

261 

254 

3 

1 

1 

1 

1 

2 

1 

1 

I 

1 

1 

1 

1 

l 

1 

1 

1 

6 

7 

8 

9 

10 

l 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

27 

34 

27 

34 

62 

-29 

-47 

21 

30 

100 

5797 

5336 

1 

4874 

5761 

4 S70 

-3894 

1 

5108 

5286 

4140 

-4103 

1 

4875 

*856 

4018 

-4713 

1 

4891 

4675 

1 

1 

1 

1 

1 

1 

1 

1 

-517 

476 

I 

1 

-504 

507 

9 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

250 

236 

219 

196 

146 

248 

235 

223 

201 

166 

247 

236 

225 

208 

191 

246 

23  7 

228 

223 

228 

246 

238 

233 

234 

247 

1 

1 

1 

1 

1 

1 

- 1 

1 

1 

1 

1 

1 

1 

1 

1 
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PEL4TIVE 

PERMEABILITY 

LOCATION  = 1 

PROBLEM 

LJ51 

J 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

29 

1 

1 

I 

1 

I 

1 

1 

1 

1 

1 

28 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

27 

1 

l 

1 

1 

I 

1 

1 

1 

1 

1 

26 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

25 

IPO 

137 

247 

943 

-4008 

-5259 

-5659 

-4046 

-4946 

5660 

24 

-137 

104 

423 

5725 

1252 

-1685 

-5774 

1720 

1720 

-5774 

23 

1 

124 

81 

1 

56 

56 

1 

52 

52 

1 

22 

i 

88 

96 

1 

55 

55 

1 

52 

52 

1 

21 

1 

88 

85 

1 

56 

56 

1 

54 

54 

1 

20 

1 

101 

100 

I 

75 

75 

1 

74 

74 

1 

1 9 

1 

1 

1 

1 

1 

1 

1 

1 

1 

l 

18 

1 

143 

66 

1 

86 

83 

1 

89 

89 

1 

1 7 

-14 

-31 

27 

516 

177 

149 

520 

220 

220 

520 

16 

60 

-70 

215 

229 

2 73 

3 70 

3 72 

387 

387 

372 

15 

154 

-22  7 

209 

244 

272 

281 

295 

301 

301 

295 

14 

224 

-211 

231 

23  7 

243 

250 

254 

256 

256 

254 

13 

215 

225 

225 

229 

232 

234 

236 

237 

23  7 

236 

12 

215 

216 

220 

223 

225 

226 

228 

228 

228 

228 

1 1 

205 

209 

212 

216 

219 

222 

224 

225 

225 

224 

10 

195 

197 

203 

209 

216 

223 

226 

228 

228 

226 

9 

148 

168 

191 

2 12 

227 

236 

241 

2,44 

244 

241 

8 

73 

160 

211 

239 

260 

275 

285 

289 

289 

285 

7 

159 

183 

249 

296 

326 

351 

3 70 

382 

382 

370 

6 

20  9 

250 

302 

3 55 

407 

433 

455 

465 

465 

455 

5 

257 

318 

391 

448 

488 

510 

517 

-517 

-517 

517 

4 

289 

376 

459 

506 

-517 

-516 

504 

510 

510 

504 

3 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 1 1 1 1 1 1 I 1 1 1 
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RELATIVE  PERMEABILITY 


LOCATION  = 1 


PROBLEM  LJ5 I 


J PI  22  23  24  25  26  27  28  29  30 

29  1 1 I I I | 1 1 | | 

28  I 1 1 1 I I I | I 1 

27  1 1 I 1 1 1 1 | | j 

261  1 1 1 1 I 1 ] |] 

25  5257  4014  945  247  137  100  62  35  27  34 

24  1685  -1252  -5724  423  104  -137  100  -30  21  -47 

23  56  57  1 81  124  1 5760  4875  1 5336 

?2  55  55  1 96  88  1 5283  5107  1 3894 

21  56  56  1 85  88  1 4855  4873  I 4103 

20  74  75  1 100  101  1 4672  4889  1 4713 

*511111Ilii| 
18  83  85  1 66  143  1 476  -517  I i 


PROBLEM  LJ5I 


RELATIVE  PEPMEAPILITY  LOCATION  = I 


J 

31 

32 

33 

34 

35 

36 

29 

1 

1 

I 

1 

1 

1 

28 

I 

I 

1 

I 

1 

1 

21 

1 

1 

1 

1 

1 

1 

2 6 

1 

1 

1 

1 

1 

1 

25 

27 

2.3 

30 

29 

24 

30 

24 

-29 

-17 

-39 

-34 

-18 

-37 

23 

5797 

1 

5619 

5745 

1 

5683 

22 

^ 66« 

1 

-3912 

-4056 

1 

3913 

2! 

4139 

1 

-3487 

-3452 

1 

3392 

20 

4018 

1 

3747 

-3570 

1 

-3530 

19 

1 

1 

1 

1 

1 

1 

18 

1 

1 

1 

1 

1 

1 

17 

1 

1 

1 

1 

1 

1 

16 

1 

1 

1 

1 

1 

1 

15 

1 

1 

1 

1 

1 

1 

14 

1 

1 

1 

1 

1 

1 

13 

1 

1 

1 

1 

1 

1 

12 

1 

1 

1 

1 

1 

1 

11 

1 

1 

1 

I 

I 

1 

10 

l 

1 

1 

1 

I 

9 

1 

1 

1 

1 

1 

1 

8 

251 

261 

269 

275 

279 

281 

7 

249 

259 

267 

2 73 

277 

279 

6 

247 

257 

265 

271 

2.75 

277 

5 

246 

255 

263 

269 

273 

275 

4 

246 

254 

261 

267 

271 

273 

3 

1 

1 

1 

1 

1 

1 

2 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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PRINT  VECTOR  POTENTIAL 

$ 

VECTOR  POTENTIALS  IN  '<B/n  LOCATION  : I PROBLEM  LJ5 1 
DIVIDE  TABLE  Bv  SCALE  '‘ACTOR 


SCALE  FACTOP  = . 

1 76E  06 

J 

1 

2 

5 

4 

5 

6 

7 

8 

9 

.10 

29 

919 

916 

905 

886 

860 

827 

788 

741 

689 

68 

28 

929 

925 

914 

896 

870 

856 

796 

749 

696 

68 

27 

986 

985 

971 

952 

924 

889 

845 

794 

756 

6? 

26 

1091 

1089 

1078 

1056 

1029 

991* 

940 

882 

815 

71 

25 

4697 

4780 

4774 

4665 

4715 

4695 

4552 

4552 

4492 

4281 

2A 

8008 

8546 

8545 

7976 

8480 

8*70 

7848 

8245 

8224 

7410 

25 

8610 

8710 

8709 

85  78 

8648 

8646 

8409 

8591 

85  8 7 

7645 

22 

8858 

8875 

8875 

8851 

8819 

8818 

8659 

8546 

8542 

7766 

21 

9045 

9040 

9040 

9025 

8995 

8989 

8870 

8711 

8687 

7876 

20 

9249 

9206 

9205 

925  8 

9175 

9154 

9125 

8922 

8786 

8012 

19 

92  55 

922  9 

92  2 6 

9242 

9199 

9175 

9150 

8969 

8807 

8015 

18 

9297 

9295 

92  92 

9295 

92  80 

92  65 

9229 

9155 

8971 

8502 

17 

9555 

9555 

9556 

95  60 

9565 

9570 

9587 

9445 

9581 

9544 

16 

9409 

9411 

9417 

9428 

9445 

94  67 

9506 

9571 

9645 

9651 

15 

9462 

9465 

94  75 

9491 

9515 

9550 

9599 

9668 

9751 

9755 

14 

9509 

9515 

952  5 

9546 

9576 

9618 

9675 

9745 

9809 

9842 

15 

9549 

9554 

9568 

9591 

9626 

9672 

9751 

9806 

9872 

9908 

12 

9580 

9586 

9601 

962  7 

9665 

9712 

9774 

9851 

9Q 1 9 

9956 

1 1 

9602 

9608 

9624 

9S51 

9689 

9759 

9805 

9880 

9948 

9985 

10 

9615 

9619 

9656 

9665 

9702 

9754 

9818 

9895 

9962 

9997 

9 

9614 

9620 

9657 

9665 

9704 

9755 

9819 

9896 

9962 

9992 

8 

8422 

8428 

8442 

8467 

8501 

8546 

8600 

8660 

8709 

8727 

7 

722  9 

7255 

7246 

7267 

7296 

7554 

75  79 

7425 

7461 

746* 

6 

6054 

6058 

6049 

6066 

6090 

6121 

6156 

6190 

6215 

6206 

5 

5072 

5074 

5079 

5087 

5098 

5112 

5126 

5138 

5142 

5124 

4 

108 

107 

107 

■ 105 

105 

101 

97 

95 

88 

82 

5 

98 

98 

97 

96 

94 

91 

8R 

85 

78 

75 

2 

95 

95 

92 

91 

89 

86 

82 

78 

75 

68 

1 

9? 

92 

91 

90 

88 

85 

81 

77 

72 
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VECTOR  POTENTIALS  IN  VB/fl 


LOCATION  : I 


PROBLEM  LJ5 1 


DIVIDE  TABLE  BY  SOALE  FACTOR 
SCALE  FACTOR  = . 1 76E  06 


J 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

29 

570 

505 

437 

367 

295 

222 

149 

75 

0 

-73 

28 

575 

509 

440 

369 

2 96 

223 

149 

75 

0 

-73 

27 

602 

529 

455 

380 

304 

228 

152 

76 

0 

-75 

26 

648 

561 

477 

393 

312 

233 

155 

77 

0 

-76 

25 

3924 

3703 

3484 

3106 

2577 

2160 

1575 

555 

1 

-552 

24 

6906 

6886 

5687 

4289 

4264 

2860 

1446 

1418 

1 

-1416 

23 

6935 

6925 

5585 

4289 

4267 

2860 

1447 

1420 

1 

-1417 

22 

6983 

6969 

5644 

4304 

4277 

2868 

1453 

1421 

1 

-1418 

21 

7058 

7017 

5677 

4341 

4284 

2879 

147? 

141  1 

1 

-1408 

20 

7212 

7014 

5739 

4463 

4222 

2896 

1568 

1325 

1 

-1323 

19 

7238 

7031 

5739 

44  70 

4226 

2896 

1569 

1325 

1 

-1323 

18 

7461 

7261 

5982 

4458 

4329 

2923 

1502 

1401 

1 

-1399 

17 

93  76 

7931 

6492 

5461 

4294 

3164 

2144 

1P67 

1 

-1065 

16 

9611 

8288 

7105 

584  7 

4652 

3488 

2316 

1158 

O 

-1  156 

15 

9750 

8524 

7260 

6041 

4821 

3608 

2404 

1201 

0 

-1199 

14 

9846 

8600 

7367 

6129 

4898 

3671 

244  5 

1222 

0 

-1221 

13 

9916 

8671 

7425 

6183 

4943 

3 705 

2469 

1234 

0 

-1233 

12 

9964 

8713 

7464 

6216 

4970 

3 726 

2483 

1241 

0 

-1240 

11 

9992 

8735 

7482 

6230 

4981 

3733 

2486 

1242 

0 

-1241 

10 

10000 

8735 

74  72 

6212 

4959 

3712 

2470 

1234 

0 

-1233 

9 

9981 

8649 

7361 

6104 

4866 

3641 

2423 

121  1 

0 

-1210 

8 

864  7 

8020 

6956 

5812 

4651 

3489 

2326 

1163 

0 

-116? 

7 

7356 

6955 

6230 

5291 

4271 

3213 

2144 

107? 

0 

-1071 

6 

6099 

5809 

5293 

4569 

3717 

2799 

1864 

929 

0 

-928 

5 

3066 

2939 

2706 

2353 

1917 

1440 

954 

477 

0 

-476 

4 

74 

66 

57 

47 

38 

29 

19 

9 

0 

-9 

3 

66 

59 

51 

43 

35 

26 

17 

9 

0 

-8 

2 

62 

55 

48 

41 

33 

25 

16 

8 

0 

-8 

1 

61 

54 

47 

40 

32 

24 

16 

8 

0 

-8 

VECTOR  POTENTIALS  IN  'i!3/M 


LOCATION  = I 


PROBLEM  LJ51 


DIVIDE  TABLE  BY  SCALE  FACTO'7 
SCALE  FACTOR  = . 1 76E  06 


J 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

29 

-147 

-22  0 

-2  93 

-365 

-435 

-504 

-569 

-631 

-688 

-740 

28 

-148 

-222 

-295 

-367 

-438 

-507 

-574 

-637 

-695 

-748 

27 

-150 

-226 

-3P2 

-3  78 

-4  54 

-528 

-600 

-670 

-735 

-793 

26 

-153 

-231 

-310 

-392 

-4  75 

-560 

-646 

-733 

-813 

-881 

25 

-1572 

-2157 

-2575 

-3104 

-3482 

-3  702 

-3923 

-4281 

-4492 

-4531 

24 

-1444 

-2354 

-4262 

-4287 

-5686 

-6885 

-6905 

-7410 

-8223 

-8244 

23 

- 1 445 

-2855 

-4265 

-4287 

-5583 

-6924 

-6934 

-7643 

-8386 

-8391 

2? 

-M50 

-2863 

-4275 

-4302 

-5642 

-6969 

-6982 

-7765 

-8541 

-8546 

21 

-1469 

-2  874 

-4282 

-4339 

-5676 

-7017 

-7057 

-7875 

-8686 

-8710 

20 

-1565 

-2893 

-4220 

-4461 

-5738 

-7013 

-721  1 

-8011 

-8785 

-8921 

19 

-1567 

-2893 

-4223 

-4468 

-5738 

-703  0 

-723  7 

-8012 

-8806 

-8968 

18 

-1499 

-2921 

-4327 

-4456 

-5981 

-7260 

-7461 

-8501 

-8970 

-9134 

17 

-2142 

-3162 

-4292 

-5460 

-6491 

-7929 

-93  74 

-9543 

-9579 

-9442 

16 

-2314 

-3486 

-4651 

-5845 

-7104 

-8286 

-9609 

-9649 

-9644 

-9570 

15 

-2402 

-3606 

-4  820 

-6039 

-7258 

-8522 

-9749 

-9754 

-9730 

-9666 

14 

-2444 

-3669 

-4896 

-6128 

-7365 

-8598 

-9845 

-9840 

-9808 

-9744 

13 

-2468 

-3  704 

-4942 

-6181 

-7423 

-8670 

-9914 

-9907 

-9871 

-9804 

12 

-2482 

-3  725 

-4969 

-6215 

-7463 

-8712 

-9962 

-9954 

-9917 

-9849 

1 1 

-2485 

-3  732 

-4979 

-6228 

-7481 

-8734 

-9991 

-9984 

-9947 

-9879 

10 

-2470 

-371 1 

-4957 

-6211 

-7471 

-8734 

-9998 

-9995 

-9961 

-9894 

9 

-2423 

-3640 

-4865 

-6102 

-7359 

-8648 

-9979 

-9991 

-9961 

-9895 

8 

-2325 

-3488 

-4650 

-5810 

-6955 

-8019 

-8646 

-8726 

-8709 

-8660 

7 

-2144 

-3212 

-42  70 

-5289 

-6229 

-6954 

-7355 

-7463 

-7460 

-7424 

6 

-1863 

-2798 

-3716 

-4568 

-5291 

-5808 

-6098 

-6205 

-6214 

-6189 

5 

-954 

-1440 

-1916 

-2352 

-2  705 

-2939 

-3065 

-3123 

-3141 

-3138 

4 

-19 

-28 

-3  8 

-47 

-57 

-66 

-7.4 

-81 

-88 

-93 

3 

-17 

-26 

-35 

-43 

-51 

-59 

-66 

-72 

-78 

-83 

2 

-16 

-25 

-33 

-40 

-48 

-55 

-62 

-68 

-73 

-78 

I 

-16 

-24 

-32 

-40 

-47 

-54 

-61 

-67 

-72 

-77 
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VECTOR  POTENTIALS  IN  ■■IB /I 


LOCATION  = I 


PROBLEM  LJ5I 


DIVIDE  TABLE  BY  SCALE  FACTOR 
SCALE  FACTOR  = .176E  06 


J 

31 

32 

33 

34 

35 

36 

29 

-787 

-826 

-860 

-885 

-904 

-915 

28 

-795 

-835 

-869 

-895 

-914 

-925 

27 

-844 

-888 

-923 

-951 

-970 

-982 

26 

-939 

-990 

-1028 

-1055 

-1077 

-1089 

25 

-4552 

-4696 

-4716 

-4665 

-4774 

-4779 

2 4 

-784  7 

-8468 

-8479 

-7976 

-8542 

-8545 

23 

-8408 

-8645 

-8647 

-8577 

-8708 

-8709 

22 

-8659 

-881  7 

-8818 

-8830 

-8874 

-8874 

21 

-8869 

-8988 

-8992 

-9022 

-9039 

-9039 

20 

-9122 

-9153 

-9172 

-9237 

-9204 

-9205 

19 

-9129 

-9174 

-9198 

-9241 

-9226 

-9228 

18 

-9228 

-9263 

-92  79 

-9292 

-9291 

-9292 

17 

-9386 

-9369 

-93  62 

-9359 

-9355 

-9352 

16 

-9505 

-9466 

-9442 

-942  7 

-94  16 

-9410 

15 

-9598 

-9549 

-9514 

-9490 

-94  74 

-9465 

14 

-9672 

-9617 

-9575 

-9545 

-9525 

-9513 

13 

-9730 

-9671 

-9625 

-9591 

-9567 

-9553 

12 

-9773 

-9711 

-9662 

-9626 

-9600 

-9585 

11 

-9802 

-9738 

-9688 

-9650 

-9623 

-9607 

10 

-9816 

-9753 

-9701 

-9663 

-9635 

-9619 

9 

-9818 

-9754 

-9703 

-9S64 

-9636 

-9620 

8 

-8599 

-8545 

-8500 

-8466 

-8442 
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Program  FLUX  and  associated  subroutines  provide  the  graphical  output  for  this 
analysis.  The  graphics  terminal  has  a hard-copy  unit  for  making  a permanent  record  of 
the  CRT  display. 

The  first  picture  (Figure  1)  provides  a record  of  the  alternator  configuration 
under  consideration  in  problem  LJ51.  Dashed  lines  represent  internal  boundaries 
which  clearly  outline  iron-air  interfaces  for  rotor  and  stator.  The  upper  and  lower 
dashed  lines  are  iron-air  interfaces  for  airspace  used  to  reduce  the  flux  through  these 
regions  to  zero.  Stator  and  damper  slots  are  shown  although  current  values  are  not 
considered  for  either  during  this  analysis. 

In  response  to  input  of  parameter  1 on  the  control  teletypewriter,  the  magnetic 
field  is  drawn  over  the  machine  outline  (Figure  2).  This  plot  is  equal  values  of  mag- 
netic induction,  B.  Series  of  closed  curves  are  drawn  which  are  symmetrical  around  the 
pole  centerline.  In  several  cates,  the  curves  close  on  themselves  across  the  vertical 
boundaries  of  the  configuration.  In  those  instances  where  curves  appear  to  fail  to 
close,  this  problem  is  really  no  indication  that  the  solution  algorithm  has  failed  but 
that  the  program  for  drawing  the  display  has  failed  at  that  point. 

In  response  to  input  of  parameter  2 on  the  control  teletypewriter,  the  magnetic 
vector  potential  is  drawn  over  the  machine  outline  (Figure  3).  The  equipotential  lines 
form  closed  curves  which  display  a symmetrical  relation  around  the  pole  centerline. 
These  curves  represent  flux  lines,  and  the  flux  lines  at  opposite  vertical  boundaries  are 
in  opposite  direction  just  as  the  sign  reversal  in  magnitudes  of  magnetic  vector  poten- 
tial maps  printed  by  the  previous  program. 

A modification  was  made  in  this  program  to  print  curves  of  small  interval  near  the 
maximum  and  minimum  values  of  the  magnetic  vector  potential.  This  modification 
tends  to  isolate  the  kernels  which  occur  in  the  field  coil  winding  on  the  rotor.  The 
record  also  shows  paths  of  leakage  flux  which  occur  in  air  between  field  poles.  Note 
the  flux  lines  at  the  iron-air  interface  are  perpendicular  as  desired  (Figures  4 and  5). 

The  number  of  lines  to  be  plotted  may  be  specified.  The  more  lines  requested, 
the  smaller  the  interval  between  equipotential  lines.  The  number  of  lines  drawn  is 
specified  and  should  not  be  interpreted  as  degree  of  magnetic  saturation. 


Figure  2.  U51  Magnetic  Field  (11  Feb  76): 
PR1NB  MOD;  PLOTB  MOD;  FLUX  MOD. 
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Figure  3.  LJ5I  Vector  Potential  (1 1 Feb  76): 
PR1NB  MOD;  PLOTB  MOD;  FLUX  MOD;  10  LINES. 
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